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Abstract 

In this work we introduce a generalization of the Jauch and Rohrlich 
quantum Stokes operators when the arrival direction from the source is 
unknown a priori. We define the generalized Stokes operators as the 
Jordan-Schwinger map of a triplet of harmonic oscillators with the Gell- 
Mann and Ne'eman SU(3) symmetry group matrices. We show that the 
elements of the Jordan-Schwinger map are the constants of motion of 
the three-dimensional isotropic harmonic oscillator. Also, we show that 
generalized Stokes Operators together with the Gell-Mann and Ne'eman 
matrices may be used to expand the polarization density matrix. By tak- 
ing the expectation value of the Stokes operators in a three-mode coherent 
state of the electromagnetic field, we obtain the corresponding generalized 
classical Stokes parameters. Finally, by means of the constants of motion 
of the classical three-dimensional isotropic harmonic oscillator we describe 
the geometric properties of the polarization ellipse 
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1 Introduction 



In both classical and quantum optics, Stokes parameters have proven to be 
intuitive and practical tools for characterizing the polarization state of light 
[1-6]. 

A classical or quantum electromagnetic wave propagates, in general, in an 
arbitrary but fixed direction in space. However for the study of the polariza- 
tion properties of the wave, the knowledge of the propagation direction of the 
wave allows us to chose a coordinate system in such away that the propaga- 
tion is along the z axis and the polarization vector rests on the x — y plane (i. 
e. it has only two polarization components) [10]. Also, the knowledge of the 
propagation direction allows to use a 2D apparatus (polarizers, wave plate ro- 
tators, etc.) placed perpendicular to the wave propagation direction to measure 
the polarization characteristics of the wave. The works [1-9] were done under 
the assumption that the arrival direction from the source of the electromag- 
netic wave was known. In [7] G. Stokes studied the polarization properties of a 
quasi-monochromatic plane wave of light in an arbitrary polarization state by 
introducing four quantities, known since then as the Stokes parameters. Wiener 
used the 2x2 unit matrix and the Pauli matrices as a basis to expand the 
coherence tensor [8]. Fano [9] showed that the Stokes parameters are the ex- 
pansion coefficients of the coherence tensor. Stokes parameters obtained under 
an a priori knowledge of the propagation direction will be refered in this work as 
the usual classical or quantum Stokes parameters, and they are well described 
in references [10, 15] and [11], respectively. 

When we do not known a priori the propagation direction of the wave, 
we already do not get the above adequate chose of the coordinate system, so 
that, in general, the three components of the polarization vector are nonzero. 
In this case, the three-dimensional coherence tensor must be used to obtain a 
complete polarization characterization [12, 13]. Roman [12] used the basis of 
nine 3x3 matrices which form the Kemmer algebra to define the generalized 
Stokes parameters as the expansion coefficients of the correlation matrix. In 
[13] Carozzi et. al. defined the generalized Stokes parameters as the expansion 
coefficients of the spectral density tensor in terms of the SU (3) Gell-Mann and 
Ne'eman matrices. 

In this work we introduce a generalization of the Jauch and Rohrlich quan- 
tum Stokes operators when the arrival direction from the source is unknown a 
priori. For simplicity we study the case of a monochromatic quantized plane 
electromagnetic wave that propagates in a fixed but arbitrary direction in space. 
Also, we will use U = u = fi = 1, where fi is the mass of each ID harmonic 
oscillator and ui is the angular frequency of either the electromagnetic wave or 
each harmonic oscillator. In section 2, we define the generalized quantum Stokes 
operators as the Jordan-Schwinger map of a triplet of harmonic oscillators with 
the Gell-Mann and Ne'eman Aj matrices of the SU(3) symmetry group. We 
show that the elements of the Jordan-Schwinger map are the constants of mo- 



2 



tion of the quantum 3D isotropic harmonic oscillator. Also, we show that the 
generalized Stokes operators together with the Ai matrices may be used to ex- 
pand the polarization matrix. In section 3, by taking the expectation value of 
the generalized quantum Stokes operators in a three-mode coherent state of the 
electromagnetic field, we obtain the corresponding generalized classical Stokes 
parameters. In section 4, by means of the classical constants of motion of the 
3D isotropic harmonic oscillator wc describe the geometrical properties of the 
polarization ellipse. Finally, in section 5, we give some concluding remarks. 

2 Jordan- Schwinger map and the harmonic os- 
cillator constants of motion 

Usual classical Stokes parameters are defined as the expansion coefficients of 
the polarization matrix [14, 15] as 

1 3 

J™ = -y^PiSi. (i) 

where Si are the four Stokes parameters, <r = 1-2x2 and o-j, i — 1,2,3, are 
the Pauli matrices. Since the <r, matrices are such that Tr{oiOj) — 25ij and 
Tr(ooOj) = 0, then 

Tri.hoa^^s,. (2) 

2.1 Usual quantum Stokes Operators 

The usual Stokes operators for a quantized plane electromagnetic wave that 
propagates along the z axis arc defined as [11] 

So = aVo<z = a\a\ + a\a,2, S\ = aVi<z = a\a,2 + a\a\, 
S2 = a) it 2a = i(— a\a,2 + a^ai), S3 = a) tr^a = a\ai — a\d2, (3) 

where aj and aj, j — 1,2, arc the creation and annihilation operators of the 
j-th harmonic oscillator defined as 

1 1 

a j = -^(xj-iPj), a j = -j=(x j +ip j ), (4) 

with [ai,aj] = [a2,<4] = 1 and 

a t = (a t, a t), (5) 

We note that equations (3) are a particular case of the Jordan-Schwinger map 
with two kinematically independent bosons [16]. 
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In the rest of the paper, the following observation is of fundamental impor- 
tance. The quantities (3) are nothing more than the constants of motion of the 
2D isotropic harmonic oscillator, with Hamiltonian H2D = a\ai + a\ a 2 + 1. In 
fact, we can show that 

[Si,H 2D ]=0, 2 = 0,1,2,3. (6) 

The commutation relations of the Stokes operators are immediately obtained 
from the properties of the Jordan-Schwinger map [16]. This leads us to the 
SU(2) Lie algebra 



St S m 
T'~2~ 



s 

i^mny, £,m,n= 1,2,3, (7) 



where ee mn is the well known totally antisymmetric tensor. 

We note that the angular momentum and the energy minus the zero point 
energy of the 2D isotropic harmonic oscillator are equal to 

L z = S2, H 2 d — 1 = Sq, (8) 

respectively. According to Jauch and Rorhlich [11], the spin of the photon is 
given by S2 and it is along the direction of propagation. Therefore, the first 
equality in (8) means that the angular momentum of the 2D isotropic harmonic 
oscillator is equal to the spin operator of the photon. 

Using equations (3), we can write the polarization matrix in terms of the 
constants of motion of the 2D isotropic harmonic oscillator (usual quantum 
Stokes operators) as 

j _\( (So) a + (S3) a (Si) a + i(S 2 )a \ , Q x 

J2D ~ 2\(S 1 ) a -i(S 2 ) a (So) a -(S 3 ) a ) W 

where (Si) a means the classical limit of the Stokes operators, which will be 
found in section 3 by taking their expectation values when the states of the 
electromagnetic field are expressed as coherent or semiclassical states. 

The physical and geometrical implications of the equality between the Stokes 
operators and the constants of motion of the 2D isotropic harmonic oscillator 
are extensively discussed in Ref. [17]. 



2.2 Generalized Stokes Operators 

When the direction of arrival from the source is unknown a priori, we generalize 
the quantum Stokes operators as follows. By using the Gell-Mann and Ne'eman 
Xi matrices of the SU (3) symmetry group [18] and the triplet of independent 
harmonic oscillators (three independent bosons) a) = (a\, a\, a\) , we define the 
generalized quantum Stokes operators as the Jordan-Schwinger map Sj = a^Xta, 
which explicitly are given by 
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S = a[ai + a\a2 + a\a3, Si = a\a2 + a\a\, S 2 = i{-a\a 2 + a\a\), 

S3 = a\a\ — a\a 2 , S4 = a\a 3 + a\a\, S5 = i(a\a\ — aja^), (10) 

S 6 = a 2 a 3 + a\a 2 , S 7 = i{a\a 2 — a\a 3 ), S 8 = —=(a\ai + a\a 2 - 2a\a-i), 

V3 

where we have used A = 13x3- 

From the commutation relations [ai,a{] = [02, o4] = [a3,ag] = 1, we show 
that the generalized quantum Stokes operators are the constants of motion of the 
3D isotropic harmonic oscillator with Hamiltonian H^d = a\a\+a\a,2+a\az + ^, 
i. e. 

\Si,H 3D ]=0, i = 0,...,8. (11) 

Also, by the properties of the Jordan-Schwingcr map [16], we show that gen- 
eralized quantum Stokes operators commutation rules satisfy the SU(3) Lie 
algebra 

= ifimn-^, £,m,n=l,...,8 (12) 

where the structure constants fe mn are totally antisymmetric under exchange 
of any two indices and are given by 

/l23 = 1, /l47 = /l56 = 

/246 = 2^ 7257 = 2> /345 = 77' (13) 

. i , _ Vs _ Vs 

7367 — — ^ 7458 — — , 7678 — —- 

A careful analysis leads us to show that the angular momentum operator 
L = f x p as well as the energy operator of the 3D isotropic harmonic oscillator 
are contained in the generalized quantum Stokes operators. Explicitly, we can 
show that 

3 

L\ = S7, L 2 = — S5, L3 = S 2 , Hsd — 2 = ^o- (14) 

Because of the first three equalities of this equation, the generalization of 
what we have written after equation (8) means that the angular momentum of 
the 3D isotropic harmonic oscillator essentially equal is to the spin operator of 
the photon. 

We generalize the definition of the polarization matrix as follows 



1 1 8 

^D = 7Ao(E ) a + -^A 1 (S,) [ „ (15) 
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again, (Ej) a is the classical limit of the generalized quantum Stokes operators, 
which as it will be shown in the next section, are the expectation values of the 
operators Ej in a coherent state of the electromagnetic field. 

Since the Ai matrices are such that Tr(XiXj) — 28{j and Tr(XoXi) = 0, 
i, j = 1, 8, then 

Tr(J 3D X 3 ) = \Tv ^AiA^Ei)^ = (£,•>„. (16) 

By using equations (10) and (15) the polarization matrix in terms of the 3D 
isotropic harmonic oscillator constants of motion (or the generalized quantum 
Stokes operators) takes the form 

^ |(E ) a + |(E 3 ) a + ^(E 8 ) Q i\{Y. 2 ) a i(£ 4 ) a - i|(E 5 ) Q 

J3D = |(El) a + *|(E 2 ) Q \ (Eo) a — \ (E 3 ) a + {^s)a \ (^e)a — i\ a 

\ a + i\{^b) a 5(S 6 ) Q + i\{Hl) a 3 (E ) a - 4j (Es)c 

(17) 

We observe that the A coefficient in equation (15) is such that equation (17) 
reduces to J 2 _d when the propagation direction of the plane electromagnetic 
wave is selected to be along the z axis. Also, we note that our definitions of 
J 3 £> is such that the trace of J 2 _d and J%d remains invariant. 

It is important to note that the polarization matrix (15) formally can be 
defined for purely quantum states. This means that it can be defined without 
taking the expectation values in a semiclassical state of the electromagnetic field 
of the Stokes operators A^. In this way, equation (16) becomes to Tr{J%r>Xj) = 
Tij. However, the implications of this definition are out of the scope of this 
work. 



3 Generalized classical Stokes parameters 

We will obtain the classical limit for the generalized quantum Stokes operators. 
To do this, we proceed as in Rcf. [20] to obtain the classical limit of the usual 
Stokes operators by taking the expectation value of the operators (3) in a two- 
mode coherent state of the electromagnetic field. In our case, we compute the 
mean value of the generalized quantum Stokes operators (10) in a three-mode 
coherent state of the electromagnetic field 

|ai,a 2 ,a 3 )= > , N,n 2 ,n 3 ). (18) 

ni,n,2, "3=0 

This leads us to the generalized classical Stokes parameters 

(E ) a = |aoi| 2 + |«02| 2 + |ao3| 2 , = 2|aoi||ao2| cos A 2 i, 
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(T,2) a = 2|aoi||ao2| sin A21, (S 3 ) a = |aoi| 2 - |ao 2 | 2 , 

(E 4 ) a = 2|aoi||ao3|cosA 3 i, (£ 5 ) a = 2|aoi||a 03 | sin A 3i , (19) 

(£ 6 ) a = 2|a 02 ||ao3| cos A 32 , (S 7 ) Q = 2|a 02 ||ao3| sin A 32 , 

(S 8 ) Q = |aoi| 2 + |ao 2 | 2 - 2|a 3| 2 , 

where a, = |aoi| exp(i0i) and Ay = fa — fa is the classical phase difference. 

As we will show in the next section, equation (19) represents the Stokes 
parameters for three classical oscillations of amplitudes |ao»| and phases fa. It 
is immediate to note that these equalities reduce to the usual classical Stokes 
parameters as a particular case when the amplitude and phase of the third 
oscillation vanishes. 



4 Classical Stokes parameters, classical 3D isotropic 
harmonic oscillator constants of motion and 
the geometric properties of the polarization 
ellipse 

The classical three-dimensional isotropic harmonic oscillator is a particle that 
moves under the force 

F = -r. (20) 

By solving the Newton's second law and imposing the initial conditions r t=0 = 
x and vt=o = v Q we obtain the solutions 

Xi = ciicost + biSint, i = 1,2,3 (21) 

where a, = x oi and bi — v oi . It is easy to see that these solutions satisfy the 
ellipsoid equation 

x\{a\ + b\ + aj + bj) + x\(a\ + b\ + a\ + fe 2 ) + x\{a\ + b\ + a\ + b\) 
-2x 1 x 2 (a 1 a 2 + b x b 2 ) - 2x 2 x 3 (a 2 a 3 + b 2 b 3 ) - 2x 1 x 3 (aia 3 + 6163) (22) 
= (0162 - a2^i) 2 + {a 2 b 3 - a 3 b 2 ) 2 + (a 3 h - aib 3 ) 2 . 

This means that the orbit of the classical 3D isotropic harmonic oscillator is 
contained in the ellipsoid. Moreover, since the classical 3D isotropic harmonic 
oscillator potential has spherical symmetry, its orbit is restricted to be in the 
perpendicular plane to the classical angular momentum L c ; = r x p. Thus, the 
elliptic orbit of the classical 3D isotropic harmonic oscillator is given by the 
intersection curve of the ellipsoid (22) and the orthogonal plane to L c /, which 
contains the origin of coordinates. 

Equation (21) can be written in oscillation form as 

Xi = \a 0i \ sm(<jjt + fa), (23) 
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with 

ai = \aoi \ sin^, b t = \a ai \ cos<&. (24) 
These equalities imply that 

a li =a 2 l + &?, sin^ = , cos<& = h% (25) 

The amplitudes and phases of the three classical oscillations of equation (23) 
depend on the initial conditions ai and bi according to equations (25). Thus, 
if we substitute equation (25) into equation (19), we incorporate the initial 
conditions in the generalized classical Stokes parameters (constants of motion 
of the classical 3D isotropic harmonic oscillator). In particular, at t = 0, the 
constant of motion of the angular momentum vector is 

L cl = «£ 7 ) a , -<E 5 )a, (£2>a) = a x b, (26) 

where a = (ai,a 2 ,a 3 ), b = (bi,b 2 ,b 3 ), and the ellipsoid equation (22) results 
to be 

4(£ ) Q - (E 8 ) Q 1 \ 2 / 4<S ) a - (S 8 ) a 1 



6 2 N ' ) £ \ 6 2 

J /2(So)a + (E 8 )a\ . . . . . . 

+X3 - XlX 2 (£l) a - X2Xz{^) a - XiXz{^i)c 



= i((S 7 )2+(E 5 )2+(E 2 )2). (27) 

Following the definition of the Euler angles of Rcf. [21], we perform a rotation 
such that the direction of the new X\ axis coincide with that of the line of nodes, 
and the direction of the new X3 axis coincide with that of L c ;. The direction 
of the line of nodes (direction of the intersection line between the orbit and the 
xi — x 2 plane) is found by a unitary vector in the x\ — x 2 plane, perpendicular 
to L ci = ((E 7 ) ce ,-(S 5 ) ce , (£ 2 ) a ). This leads us to 

• a ^ a 
smcp = n x = ±- 



y/(E 7 )l + (E 5 )| 
/■£ ) 

cos 4> = n y = T ; ■ (28) 

The orthogonality between L c / and the ellipse plane leads to 

cos 6 = (S2) " . (29) 

^((S 7 )2 +(E 5 )2 +(S 2 )2) 

On the other hand, it is well known that the constants of motion of the classical 
3D isotropic harmonic oscillator, in addition to the energy and the angular 
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momentum are given by the symmetric Runge-type tensor [22] 

A ij ; = ^(PiPj +uj 2 x l x j ) 7 i,j = 1,2,3 (30) 

It can be shown that the contraction of this equation with the components of 
L c ; yields to zero. This means that all the geometric characteristics of the orbit 
must be determined by A^. In fact, we can show that 

, 2(S ) a + (£ 8 ) a (S 3 )a . 2(S ) a + (£ 8 ) a 

2^11 = 7. 1 7, ) 2A 22 = 7. 

6 2 6 

(So) a - (S 8 ) a 1 

<^33 — g , ZA 12 — 1 ' a ' 13 " 

2A 23 = l(E 6 ) a . 

This shows that the geometric properties of the polarization ellipse exclusively 
depend on the generalized classical Stokes parameters. It can be shown that the 
eigenvalues of A^j depend only on the energy and the magnitude of the angular 
momentum of the classical 3D isotropic harmonic oscillator [22]. A^ has an 
eigenvector in the direction of the angular momentum, and the other two eigen- 
vectors are in the directions of the minor and major axis of the elliptical orbit 
[22]. Also, the eigenvectors of a symmetric rank two tensor are determined by 
its eigenvalues as well as its components [23] . All this leads us to conclude that 
the principal axis directions on the elliptical orbit are completely determined 
by the generalized classical Stokes parameters. Also, since L c / is orthogonal to 
the polarization ellipse, then L c ; points along the propagation direction of the 
electromagnetic wave. 



(S3)a 
2 ' 

^<£4>a, (31) 



5 Concluding Remarks 

This work links the quantum optics to classical optics by means of quantum me- 
chanics and it is a useful extension of the generalized classical Stokes parameters 
into the quantum domain. 

Although already there are treatments of the classical Stokes parameters in 
the case of an a priori unknown of the electromagnetic propagation [12, 13], 
our treatment results to be novel in the following aspects. We have introduced 
a generalization of the quantum Stokes parameters of Jauch et. al. [11] using 
the Jordan-Schwinger map, three independent bosons and the Gell-Mann and 
Ne'eman SU(3) symmetry group matrices. It was shown that the generalized 
quantum Stokes operators result to be the expansion coefficients of the polar- 
ization matrix in terms of the Gell-Mann and Ne'eman SU (3) matrices. The 
semiclassical limit of the generalized Stokes operators were achieved by taking 
their expectation values in a three-mode coherent state of the electromagnetic 
field. We have shown that the resulting generalized classical Stokes parameters 
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are consistent with the generalized classical Stokes parameters which have re- 
cently been published [13]. Thus, our treatment is more general than those of 
references [12, 13], which are restricted to the classical aspects of electromag- 
netic polarization. 

We described by means of the classical 3D isotropic harmonic oscillator con- 
stants of motion the geometrical properties of the polarization ellipse. Partic- 
ularly, we showed that the ellipsoid coefficients and the symmetric Runge-type 
tensor of the classical 3D isotropic harmonic oscillator are completely deter- 
mined by the generalized classical Stokes parameters. Also, we showed that the 
first two Euler angles are intimately related to the components of the orbital 
angular momentum of the classical 3D isotropic harmonic oscillator. 

Finally, we emphasize that the number of independent generalized classical 
Stokes parameters are six. This is because in going from (22) to (27), all of 
them were written in terms of the six parameters, <n and &i, i = 1,2,3 which 
contain the initial conditions of the classical 3D isotropic harmonic oscillator. 
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